Search for: All records

Abstract Fingering convection (also known as thermohaline convection) is a process that drives the vertical transport of chemical elements in regions of stellar radiative zones where the mean molecular weight increases with radius. Recently, Harrington & Garaud used three-dimensional direct numerical simulations (DNS) to show that a vertical magnetic field can dramatically enhance the rate of chemical mixing by fingering convection. Furthermore, they proposed a so-called “parasitic saturation” theory to model this process. Here, we test their model over a broad range of parameter space, using a suite of DNS of magnetized fingering convection, varying the magnetic Prandtl number, magnetic field strength, and composition gradient. We find that the rate of chemical mixing measured in the simulations is not always predicted accurately by their existing model, in particular when the magnetic diffusivity is large. We then present an extension of the Harrington & Garaud model which resolves this issue. When applied to stellar parameters, it recovers the results of Harrington & Garaud except in the limit where fingering convection becomes marginally stable, where the new model is preferred. We discuss the implications of our findings for stellar structure and evolution. 

ABSTRACT Zahn’s widely used model for turbulent mixing induced by rotational shear has recently been validated (with some caveats) in non-rotating shear flows. It is not clear, however, whether his model remains valid in the presence of rotation, even though this was its original purpose. Furthermore, new instabilities arise in rotating fluids, such as the Goldreich–Schubert–Fricke (GSF) instability. Which instability dominates when more than one can be excited, and how they influence each other, were open questions that this paper answers. To do so, we use direct numerical simulations of diffusive stratified shear flows in a rotating triply periodic Cartesian domain located at the equator of a star. We find that either the GSF instability or the shear instability tends to take over the other in controlling the system, suggesting that stellar evolution models only need to have a mixing prescription for each individual instability, together with a criterion to determine which one dominates. However, we also find that it is not always easy to predict which instability ‘wins’ for given input parameters, because the diffusive shear instability is subcritical, and only takes place if there is a finite-amplitude turbulence ‘primer’ to seed it. Interestingly, we find that the GSF instability can in some cases play the role of this primer, thereby providing a pathway to excite the subcritical shear instability. This can also drive relaxation oscillations, which may be observable. We conclude by proposing a new model for mixing in the equatorial regions of stellar radiative zones due to differential rotation. 

The solar tachocline is a shear layer located at the base of the solar convection zone. The horizontal shear in the tachocline is likely turbulent, and it is often assumed that this turbulence would be strongly anisotropic as a result of the local stratification. What role this turbulence plays in the tachocline dynamics, however, remains to be determined. In particular, it is not clear whether it would result in a turbulent eddy diffusivity, or anti-diffusivity, or something else entirely. In this paper, we present the first direct numerical simulations of turbulence in horizontal shear flows at low Prandtl number, in an idealized model that ignores rotation and magnetic fields. We find that several regimes exist, depending on the relative importance of the stratification, viscosity and thermal diffusivity. Our results suggest that the tachocline is in the stratified turbulence regime, which has very specific properties controlled by a balance between buoyancy, inertia, and thermal diffusion.